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ABSTRACT 

Sequence analysis of the 3' part (8 kb) of the 
polymerase gene of the torovirus prototype Berne virus 
(BEV) revealed that this area contains at least two open 
reading frames (provisionally designated ORFIa and 
ORFIb) which overlap by 12 nucleotides. The complete 
sequence of ORFIb (6873 nucleotides) was 
determined. Like the coronaviruses, BEV was shown 
to express its ORFIb by ribosomal frameshifting during 
translation of the genomic RN A. The predicted tertiary 
RNA structure (a pseudoknot) in the toro- and 
coronaviral frameshift-directing region is similar. 
Analysis of the amino acid sequence of the predicted 
BEV ORFIb translation product revealed homology 
with the ORFIb product of coronaviruses. Four 
conserved domains were identified: the putative 
polymerase domain, an area containing conserved 
cysteine and histidine residues, a putative helicase 
motif, and a domain which seems to be unique for toro- 
and coronaviruses. The data on the 3' part of the 
polymerase gene of BEV supplement previously 
observed similarities between toro- and coronaviruses 
at the level of genome organization and expression. 
The two virus families are more closely related to each 
other than to other families of positive-stranded RNA 
viruses. 

INTRODUCTION 

Berne virus (BEV) is the prototype of the proposed family 
Toroviridae, a recently characterized group of enveloped, 
positive-stranded RNA viruses (1) which have been found in 
association with enteric infections in horses, cattle and man (2). 

The BEV genome consists of a single, infectious, 
polyadenylated RNA molecule (3). From U.V. transcription 


mapping data its size was calculated to be between 25 and 30 
kilobases (kb; 4). The genome organization and expression of 
BEV are comparable to those of the coronaviruses (for a review 
see 5). Sequence analysis has revealed at least five open reading 
frames (ORFs) on the BEV genome (4). These ORFs are 
expressed from a nested set of mRNAs, each containing one of 
the ORFs at their 5' end. The BEV mRNAs are 3' coterminal, 
but, in contrast to the coronaviral mRNAs, no common leader 
sequence has been detected at their 5' end (4). 

In addition to similarities at the level of genome organization 
and expression, the presence of amino acid sequence homologies 
in replicase proteins is considered indicative for ancestral 
relationships between groups of viruses (6, 7). Recently, we have 
demonstrated that the large ORF which is located upstream of 
the BEV peplomer protein gene is only present on genomic RNA. 
Therefore this ORF most likely encodes the viral RNA-dependent 
RNA polymerase (POL) or replicase (4). Since the termination 
codon of this ORF maps at 6603 nucleotides (nt) upstream of 
the poly(A) tail, the complete replicase encoding region of BEV 
is extremely large (18 to 23 kb). Among positive-stranded RNA 
viruses, its size is comparable only to the POL gene of the 
coronaviruses (8, 9). 

The coronaviral POL gene has been shown to contain two 
ORFs (8) of about 12 kb (ORFIa, formerly named FI) and 8 
kb (ORFIb, formerly named F2). The 3' end of ORFIa overlaps 
with the 5' end of ORFIb, which is expressed by ribosomal 
frameshifting. This process involves RNA sequence elements in 
the ORFla/ORFlb overlap region, including an RNA pseudoknot 
( 10 , 11 , 12 ). 

In this paper the sequence analysis of a substantial part of the 
BEV polymerase gene is reported. Both the primary structure 
of the BEV polymerase and the mechanism used to express the 
3' part of the POL gene indicate that toro- and coronaviruses 
are indeed evolutionarily related. 


To whom correspondence should be addressed 

Present address: Institute of Medical Microbiology, Department of Virology, State University of Leiden, Postbus 320, 2300 AH Leiden, The Netherlands 


Downloaded from http://nar.oxfordjournals.org/ at University of Exeter on July 22, 2015 







4536 Nucleic Acids Research, Vol. 18, No. 15 


MATERIALS & METHODS 

cDNA synthesis and cloning 

The preparation of a randomly primed cDNA library from 
intracellular poly(A)-selected BEV RNA and the selection of 
BEV-specific cDNA clones have been described previously (4). 

Sequence analysis and computer analysis of sequence data 

Nucleotide sequence analysis and assembly of sequence data were 
performed as described previously (4). Amino acid sequence 
similarity searches were carried out using the FASTA program 
(13) and the NBRF protein identification resource (release 20.0). 
Dot matrix comparisons, sequence alignments and RNA 
secondary structure predictions were generated using the 
COMPARE, GAP and FOLD options from the software provided 
by the Computer Genetics Group/University of Wisconsin 
(version 5, 1989; 14). 

Construction of pBSFS 

To study ribosomal frameshifting in the ORFla/ORFlb 
overlapping region the expression vector pBSFS was constructed. 
The unique EcoRI site of plasmid pBS (Stratagene) was digested 
and sticky ends were filled using the Klenow fragment of DNA 
polymerase I (15). An Xhol linker (5' CCTCGAGG 3') was 
inserted into the filled EcoRI site, resulting in pBSXho. An Xhol 
DNA fragment containing the matrix protein (M) gene of mouse 
hepatitis virus strain A59 (MHV-A59) (a gift from Peter Rottier) 
was cloned into the Xhol site of pBSXho to generate construct 
pBSM. pBSM was digested with Smal and Sail, sticky ends were 
filled as described above, and the vector was religated. In this 
way the BamHI site in the multiple cloning region of pBSM was 
deleted. Subsequently, Bglll (5' CAGATCTG 3') and BamHI 
(5' CGCGGATCCGCG 3') linkers were inserted into the filled 
AccI and Styl sites, respectively, of the MHV-A59 M gene in 
pBSM (nt positions 656 and 693 in the sequence as presented 
in Fig. 1 of reference 16), which resulted in construct pBSMBB. 
The cDNA insert of BEV clone 129 (Fig. 1) was digested with 
Sau3A and a 675 bp fragment (nt positions 563 to 1237 in Fig. 
2) was purified. pBSFS was generated by cloning this Sau3A 
fragment, which contains the ORFla/ORFlb overlapping region, 
between the digested and dephosphorylated BamHI and Bglll sites 
of pBSMBB. The orientation of the insert and the nucleotide 
sequences at the MHV M/BEV POL junctions were examined 
by sequence analysis. The hybrid gene in pBSFS is under the 
control of the pBS T7 promoter. 

In vitro transcription and translation 

Plasmid DNA of construct pBSFS was linearized by digestion 
with HindHI. In vitro transcription, using T7 RNA polymerase, 
and in vitro translation of the resulting RNA transcripts in a rabbit 
reticulocyte lysate were performed as described previously (17). 

In vivo expression of pBSFS 

Hela cells (lxlO 6 ) were infected (m.o.i. of 5 p.f.u.) with 
recombinant vaccinia virus vTF7 —3, which contains the T7 
polymerase gene under the control of a vaccinia promoter (18). 
After a 90 minute incubation, the cells were transfected with 
pBSFS (5 fig) as described by Gorman (19). The culture 
supernatant was removed at 13 Vi hours post infection (p.i.) and 
replaced by methionine-free medium. A 30 minute labelling with 
60 /zCi/ml [ 35 S]methionine was performed at 14 hours p.i. Cells 
were lysed in 300 fi\ lysis buffer (10 mM Tris-HCl pH 7.4, 1 
mM EDTA, 150 mM NaCl, 1% Triton X-100). Lysates were 


clarified by centrifugation at 10,000 g for 60 minutes (4°C) and 
stored at — 70°C. 

Immunoprecipitation and analysis of pBSFS translation 
products 

To immunoprecipitate pBSFS translation products, a monoclonal 
antibody (J. 1.3; a gift from John Fleming and Stephen Stohlman; 
20) directed against the N-terminal region of the MHV-A59 M 
protein and an antipeptide serum raised against the C-terminal 
18 amino acids (aa) of the same protein (Rottier, manuscript in 
preparation) were used. Antiserum and SDS (see below) were 
added to in vitro translation mixtures or transfected cell lysates 
and the mixture was diluted in lysis buffer to give a final volume 
of 300 fi\. Following an overnight incubation at 4°C, 50 u 1 of 
a 3M KC1 solution and 25 /d Pansorbin (Calbiochem) were added. 
After an additional hour at 4°C, immune complexes were spun 
down and washed three times in lysis buffer (500 /zl). Throughout 
the entire procedure SDS was present at a concentration of 0.1 % 
(N-terminal serum) or 1 % (C-terminal serum). Immuno- 
precipitated proteins were suspended in Laemmli sample buffer 
and analyzed on 12.5% SDS-polyacrylamide gels. 

RESULTS 

Partial cloning and nucleotide sequence analysis of the BEV 
POL gene 

The cloning and sequence analysis of about 10 kb of the BEV 
genome (starting at the poly(A) tail) and the identification of about 
3 kb of the putative BEV POL gene have been reported previously 
(4). To extend the BEV sequence towards the 5' end of the viral 
RNA, cDNA clones were selected from a random-primed cDNA 
library. Clones covering an additional 5 kb of sequence were 
obtained (Fig. 1) and used to determine the sequence of a 
substantial part of the BEV polymerase gene (Fig. 1 and 2). 

Analysis of the obtained POL sequence (with a total length 
of about 8 kb) revealed that the 3' part of the BEV POL gene 
contains two ORFs. At this moment we cannot exclude the 
presence of additional ORFs in the 5' part of the gene. However 
to facilitate the comparison with the coronaviral POL gene (see 
below), the presently described ORFs have provisionally been 
designated ORFla and ORFlb (Fig. 1 and 2). 

The 6873 nt ORFlb, which potentially encodes a 261K protein 
(2291 aa), covers the 3' part of the POL gene (Fig. 1 and 2). 
The first potential translation initiation codon in ORFlb is located 
at 174— 176 nt downstream of its 5' end. The 3' terminal 1053 
nt of ORFla were determined (Fig. 2) and this ORF was found 
to overlap with ORFlb over a distance of 12 nt, the latter being 
in a —1 reading frame with respect to ORFla. 

Analysis of the deduced amino acid sequence of the BEV 
ORFlb product 

A dot matrix comparison between the deduced amino acid 
sequences of the BEV ORFlb product and the corresponding 
sequences of the coronaviruses infectious bronchitis virus (IBV, 
strain Beaudette (M42); 8) and mouse hepatitis virus (MHV, 
strain A59; 12) revealed a number of domains with significant 
sequence similarity (Fig. 3). The overall identity between the 
products of ORFlb of BEV and either IBV or MHV is 
approximately 19%, which is much lower than the 56% observed 
when IBV and MHV are compared (12). However, four domains 
seem to be well conserved among toro- and coronaviruses and 
contain a much higher percentage of identical amino acids (see 
Fig. 4 and Discussion). 


Downloaded from http://nar.oxfordjournals.org/ at University of Exeter on July 22, 2015 


Nucleic Acids Research, Vol. 18, No. 15 4537 


131-— ----- -123 

129 - 120 

128 - - - - - 118 

126- - 

012345678kb 

_i_i_ J _ l _ I _i_i_ 1 _ i _ 

VPVBP BB W 


13466 E E H H H N 


ORF la 

66 

03 


ORF 1b 



P 


Fig. 1. Cloning and sequencing of the BEV ORFlb region. A number of restriction 
sites used in sequence analysis is indicated: B = BglU, E = EcoRI, H = Hindin, 
N = Ndel, P = PstI, and V = PvuII. The numbered lines indicate the position 
of the most important cDNA clones; dashed lines indicate the sequenced areas. 
Open boxes represent open reading frames. The ORF marked P is the BEV 
peplomer gene (4). Numbers indicate the distance from the start and end of ORFlb 
to the start of the poly(A)-tail. 


Detailed analysis of these well-conserved sequence motifs 
revealed that two of them (domains 1 and 3 in Fig. 4) are common 
to RNA polymerases. Domain 1 represents the so-called ‘GDD’ 
motif (21). It has been shown to be present in viral RNA- 
dependent RNA polymerases and RNA-dependent DNA 
polymerases encoded by retroid elements (22) and is thought to 
be a ‘polymerase module’. As in the coronaviral polymerase 
motif, a serine residue precedes the two conserved aspartic acid 
residues in the core of the BEV domain. The percentage of 
identical amino acids in domain 1 of BEV and IBV/MHV is 
45 — 50% and includes a 13 residue perfect match 
(KPGGTSSGDATTA) at about 70 aa upstream of the SDD core 
of the motif. 

The amino acid sequence indicated as domain 3 (Fig. 4) is 
known as the ‘NTP-binding’ or ‘helicase’ motif. Putative NTP- 
binding domains have been identified in proteins of both viral 
and cellular origin (23, 24). The viral NTP-binding proteins are 
all thought to be involved in genome replication and/or 
transcription (possibly by exerting an NTP-dependent RNA or 
DNA duplex unwinding activity). The putative BEV helicase 
domain (Fig. 4) is most similar to the coronaviral one. It is 
remarkable that the order of the polymerase and helicase domains 
in both toro- and coronaviral polymerases is reversed when 
compared to the polymerases of other positive-stranded RNA 
viruses. 

The ORFlb domain indicated as motif 2 (Fig. 4) contains a 
striking number of conserved cysteine and histidine residues. On 
the basis of the IBV sequence, the formation of three metal¬ 
binding ‘finger’ structures (25) in this region was proposed (26). 
However, comparison of the ORFlb sequences of BEV, MHV 
and IBV indicates that the metal-binding cysteine and histidine 
residues are conserved for only one of the three proposed finger 
structures. Consequently, the formation of only one conserved 
metal-binding finger in domain 2 remains possible, implying that 
the proposed ‘multi-finger’ model for IBV (26) requires 
modification. 

The fourth conserved domain is located near the C-terminus 
of the ORFlb product (Fig. 4) and seems to be unique for toro- 
and coronaviruses since no significant amino acid sequence 
similarity with other proteins could be detected. 


Analysis of the ORFla/ORFlb overlapping region 

In addition to the amino acid sequence homologies discussed 
above, our data indicated that the mechanism of expression of 
BEV ORFlb could also be comparable to the coronaviral system: 
the 3' part of the POL gene had been found to contain two 
partially overlapping ORFs and no BEV subgenomic mRNA of 
about 13.5 kb, which would be required for the independent 
translation of ORFlb, had been detected in BEV-infected cells 
(3, 4). 

Recent studies on the gene expression of retro- (27) and 
coronaviruses (10, 11, 12) have identified sequence elements and 
RNA structures which can promote ribosomal frameshifting. 
First, certain heptanucleotide RNA sequences, the so-called 
‘slippery’ sequences, have been shown to function as the actual 
frameshifting site. The heptanucleotide sequence 5' U UUA AAC 
3' has been identified as the ‘slippery’ sequence involved in IBV 
ORFla/ORFlb ribosomal frameshifting (11); it is also present 
in the corresponding region of the MHV genome (12). The same 
sequence was found from 10 to 4 nt upstream of the ORF la 
termination codon in BEV (Fig. 2, nt positions 1045—1051, and 
Fig. 5). 

Second, in the case of IBV (11) it has been shown that the 
formation of an RNA pseudoknot (28) immediately downstream 
of the ‘slippery’ sequence, is required for efficient frameshifting. 
Similar tertiary RNA structures are predicted for the frameshift- 
directing regions of MHV (12) and a number of retroviruses (11). 
A potential stem-loop structure was also identified just 
downstream of the BEV ‘slippery’ sequence (Fig. 5) and two 
sequence elements which can participate in the formation of an 
RNA pseudoknot were detected: nucleotides number 1099— 1103 
(pseudoknot 1; PK1) and 1157-1161 (pseudoknot 2; PK2) form 
stretches of 5 nt which are complementary to a part of the loop 
of the BEV hairpin (Fig. 2 and 5). From the data presented above, 
it was clear that the BEV ORFla/ORFlb overlapping region 
contains all elements required for efficient ribosomal 
frameshifting. 

The BEV ORFla/ORFlb overlapping region directs 
ribosomal frameshifting in vitro and in vivo 

To provide experimental evidence that the BEV ORFla/ORFlb 
overlapping region directs ribosomal frameshifting, construct 
pBSFS was generated. In pBSFS a restriction fragment from BEV 
clone 129 (Fig. 1) containing the ORFla/ORFlb overlap was 
cloned into a copy of the MHV-A59 26K membrane protein (M) 
gene, which is under the control of the T7 promoter (Fig. 6A). 
The 5' end of the MHV M ORF in pBSFS is fused in frame 
to the 3' end of BEV ORF la. Termination of translation at the 
ORF la termination codon will give rise to a 41K fusion protein 
(Fig. 6A). Since the 5' end of ORFlb was cloned in frame with 
the 3' end of the MHV M ORF, ribosomal frameshifting should 
result in the synthesis of an additional fusion protein of 50K. 
The two predicted products can be identified by 
immunoprecipitation with antibodies directed against the N- and 
C-terminal part of the MHV M.protein (Fig. 6A). 

Ribosomal frameshifting was studied in vitro by translation of 
pBSFS transcripts in a rabbit reticulocyte lysate. An in vivo 
analysis was performed by transfection of vaccinia virus 
recombinant vTF7—3-infected Hela cells with pBSFS DNA. The 
results of both experiments were identical: translation of pBSFS 
RNA indeed led to the synthesis of proteins of 41K and 50K. 
Both products were recognized by the antiserum directed against 
the N-terminus of the MHV M protein, whereas only the 50K 
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1 VRYITLNLTOHLQKQSIKPRRRIQYRQQGVRLGGVNLGEHQAFSNEIISTVGYTTWVSST 

1 AGTACGTTATAnACAnAAATlTAACAGACCATTTACAAAAACAAAGCCTTAAACCACGTAGACGACtTCAGtATAGACAACAGGGTGTTAGATTGGGAGGTGTTAATTTAGGrGAACATCAAGCCTTTTCAAATGAGTTAATTTCAACTGTtGGTTATACTACATGGGTTtCATCAAC 

6i vcr'dnt'hkhphevoipvnekdpe'heh'hhtolkd'nqhvvdlkpt'hhl'vradtge'qlf'alsl 

181 TGTGTGTAGAGATAATACTCACAAGCATCCATGGTTTGTTCAGATACCTGTTAATGAAAAAGATCCTGAATGGTTTATGCACAACACTCAGCTGAAAGATAATCAGTGGGTAGTTGACCTCAAACCAACTCACTGGTTGGTTAATGCTGATACTGGTGAACAACTTnTGCATTGTCCTT 

121 TDEQALKAEAILQKWSP ITODVECKFKDLKGYYTVSGFOPLWP VCp'vHiCNVR LDPVFKP 

361 AACTGATGAACAGGCTTTGAAAGCAGAGGCCATTTTGCAGAAGTGGTCACCAATAACACAGGATGTTGAGTGTTGGTTTAAAGATTTAAAAGGCTATTATACAGTTAGTGGATTTCAACCACTGTGGCCTGTTTGTCCTGTTAATATTTGTAATGTTAGACTTGATCCAGTTTTTAAGCC 

181 Q S i' V Y A D D P T H F l’ S L P V V N K H F L - A A F - Y D L Q E G f' P G K - K Q V A P H i' S L T - H L K L S D e’ D I E \ V E D 

541 TCAGTCTATAGTTTATGCTGATGATCCTACTCArmCTTTCTTTACCTGTTGTAAATAAAAATTTrCTTGCTGCATTTTATGATTTACAGGAGGGTTTTCCAGGTAAGAAACAAGTTGCACCACATATTTCTCTCACAATGTTGAAGTTGAGTGATGAAGATATAGAGAAAGTTGAAGA 

241 ild’ehv'lpnshvt’itnphhhgkh'yvc'dveglds'lhd'evvsvlr'ehg'iacdqrr’iwk'phlt 

721 CATTCTTGATGAGATGGTTTTGCCTAATTCATGGGTAACCATAACAAACCCTCATATGATGGGTAAACATTATGTTTGTGLATGTTGAGGGTCTTGATTCACTTCATGATGAGGnGTGTCGGTGCTTAGAGAGCATGGCATTGCTTGTGACCAAAAACGCCTTTGGAAACCCCATTTGAC 

301 ige' lkd'vsfokfk’dfa’isckleocdf’vklgapk'aha'ryefitt'lplgolnc* 

1 *FKLlRGAIfSACR 

901 tattggtgagttaaatgatgtaagttttgataaatttaaagattttgcaattagttgtaagttagaggattgtgattttgttaaattgggtgcgcctaaggcaaatgctaggtatgaatttataacaactcttcctttgggtgatttaaactgttgagaggtgcctggagcgcctgcagg 

13 HLCF0NGAY9SSRSKHYIDLATEYNAGIVKVNKSNTHSVE Y9S KRFMIKR VKD9SEFALA 
1081 CATCTCTGTTTTCAAAATGGCGCA1ACCAGTCTTCAAGGTCAAAACATTATATTGATTTGGCAACTGAGTATAATGCAGGCATTGTGAAAGTTAATAAATCTAATACTCATTCAGTTGAATATCAGAGTAAACGCTTTATGATTAAGCGTGTGAAAGATCAAAGTGAATTTGCATTAGCA 

73 KTAFLP5IIPHHHEK8NGEw'fEIRGPTS9WSLGOLVYAIRLGD9DYLSECGFVFNPSRDe’ 
1261 AAGACCGCGTTTTTACCATCAATAATACCTCACCATATGGAAAAACAAAATGGTGAATGGTTrTTGATACGCGGTCCTACATCTCAGTGGAGTCTTGGT&ACCTTGTCTATGCAATTTGGCTTGGTGATCAGGATTATTTGTCTGAGTGTGGTmGTTTTTAATCCATCAAGAGATGAG 

133 FLDDANORSFLANLLEPAIL KFs'HrYifqVKHCKVPYKLTL’DRVDLNGgLYDFGOYPCPNS 

1441 TTTTTGGATGATGCTAATCAAAGGTCATTTTTGGCAAATTTGCTTGAGCCTGCAATTTTAAATTTTTCTCATATCTATTGGCAAGTTAAAATGTGTAAAGTTCCCTACAAGTTAACATTGGATAATGTGGATTTGAATGGTCAGCTTTATGATTTTGGTGATTATCCTTGTCCAAATAGT 

193 vdn9Salfvlaevwshtrrp'fpvafarllaneheiptdy9hff’9nillsg'sylokalcln' 

1621 gtagataatcaatctgcgctttttgttttggctgaagtttggtctatgactcgcagaccatttccagtagcttttgcaaggcttcttgctaatgagatggaaataccaacagactatcagatgttttttcagaacattttattgtctggtagctacttggataaagcactttgtttgaat 

253 n v r p f l sdpa»ltttpffsq'hng'vwthfyn’piy’glvecrl’oef’aelp£vl'qql'vtv9gpi' 

1801 AATGTAAGACCGTTTTTGAGCGATCCAGCCAATTTAACAACTACACCATTTTTTTCACAACATAATGGTGTTTGGACACATTTTTATAACCCCATTTATGGTTTGGTTGAATGTAAFTTAGATGAArrTGCTGAATTGCCTGAAGTTTTACAACAACTTGTAACTGTGCAGGGCCCAATA 

313 TNHMTPAISVGEGVYAANVPSASATKQKIPFYDVGLC9EL TDAGVDCGEAFKYFYYLSNP 

1981 ACTAATAATATGACACCGGCTATATCAGTTGGAGAAGGTGTCTATGCAGCTAATGTTCCATCCGCATCTGCTACTAAACAGAAGATACCTTTTTATGATGTTGGTnGTGCCAAGAATTGACTGATGCTGGTGTTGATTGTGGTGAAGCATTTAAATATTTTTATTATTTATCTAACCCA 

373 AGALADVCYYOYq’GrGFYSp’KLL'AGVYDFHKRVTECYRIN’ERF TYEgAKp'RKS'SNGiHIT' 

2161 GCTGGCGCATTGGCTGATGmGCTAlTATGATTArCAGGGTACTGGGTTTTATTCCCCAAAGTTATTAGCAGGTGTTTATGATTTTATGAAAAGGGrAACAGAATGTTATAGGATAAATGAGAGATTTACATATGAGCAAGCAAAGCCAAGAAAGTCATCAATGGGAATAAATATAACA 

433 GYOQDAVYRALGPEHIARLFEYAqKAPLPF CTKIITKFALSAKARARTVSSCSFlASTIF 

2341 GGATATCAACAAGATGCAGTTTATCGTGCTTTGGGCCCTGAAAATATTGCAAGGTTGTTTGAATATGCACAGAAAGCTCCATTGCCATTTTGTACAAAAATAATAACAAAGTTTGCArrGTCAGCTAAAGCAAGAGCAAGAACAGTTTCCTCATGCTCATTTATAGCATCAACAATATTT 

493 rfahkpvtskmve'vaonsgg’fcl'igvskyg’lkf'skflkdk’yga’iegfovf’gsd’ytkcdrt’ 
2521 agatttgcacataaaccagtaacatctaaaatggttgaagtagcacagaattcaggtggcttttgcctcattggtgtatccaagtatggtttaaagttttcaaagtttttgaaagataagtacggtgctattgagggttttgatgtgtttggtagtgattacaccaagtgtgaccgtaca 

553 fplsfraitaallyelgewd’ekswiylnevnsyalotmlcdgnllnkpgg'tssgdattah 

2701 TTTCCACrGrCATTTAGAGCTTTGACAGCCGCACTTCTCTATGAATTAGGGGAGTGGGATGAGAAGTCTTGGTTATATCTTAATGAGGTGAATTCCTACATGTTGGATACTATGCTGTGTGATGGAATGCTACTTAACAAACCTGGAGGGACATCATCTGGTGATGCCACAACTGCCCAC 

613 snt‘fyiiyhvh‘yvv'afktils‘dls‘egnkvmr'iaa'hnayttg’dy9'vfrtlle'o9f’9Tnyfln' 
2881 TCAAATACTTTTTATAAnATATGGTGCATTATGTTGTAGCTTTTAAAACAArrTTGTCAGATCTTTCAGAAGGCAATAAGGTTATGAGAATTGCAGCACATAACGCCTATACAACTGGCGAnATCAAGTGTTTAATACGCTGCTTGAAGATCAAnTCAGACAAATTATTTTTTAAAT 

673 fl s'd d s f i f s' k p e ’a l k i f t c’ e n f’s n k l q t i' l h t ’k v d q t k s’ r s t 'k g h i e e f' c s a ’h i i k t d g' 
3061 tttttgtctgatgattcttttatattttcaaaacctgaggccttaaaaatttttacttgtgaaaatttttcaaataagctacaaactattttacatacaaaagtagaccaaacaaaatcatggtctacgaaaggtcatattgaggaattttgttcagcacatataattaaaacagatgga 

733 E Y H F L PSRGRLLA'sLLiLDK'LSOVDIYYHR'FVAIlCESAV'YSR YOPEFFN'GLF’oVFLORv' 
3241 GAGTATCACTTTTTACCATCTAGGGGCAlGATtGTTGGCrTCATTGCTGATTCTTGATAAGCTTTCAGATGTTGATAmATTACATGAGATTTGTTGCAATCTTGTGTGAATCTGCrGTTTATTCAAGATATCAGCCTGAGTTTTTTAATGGTTTGrnCAAGTGmCTTGATAAAGTA 

793 09f"rkdyccd'pcpp9Llere'fye’nlvftsii'sev'givdcyl'enf’kl9Cefk'9Qanfokvcf 
3421 CAACAGTTTAGAAAGGATTATTGCTGTGATCCGTGTCCACCTCAACTClTAGAGCGTGAGTTCTATGAGAATTTAGTTTTTACATCAAATAGTGAGGrGGGTATTGTTGATTGCTATCTAGAGAATTTTAAGCTTCAATGTGAGTTTAAACAACAGGCAAATTTTGATAAAGTGTGTTTT 

853 ccpnpavsvceec yvplplc"aycyyvhvvi"snhskvedkf"kcf’cg9Dnirelyivlnnsi 

3601 tgttgtcccaatcctgcagtttcagtttgtgaagagtgttatgttccccttcccctttgtgcgtattgttattatgtgcatgttgtaattagtaatcatagtaaagtagaggataagtttaaatgtttttgtggtcaagataatattagggagttgtatatagtattaaataatagtatt 

913 chy’9Ckhcvesdr'lrislls'ovo'9ivrlpg'fksnsasiak'hgv'a9lltsv'okvdvsldwn 

3781 TGTATGTATCAGTGTAAAAATTGTGTAGAGAGTGATAGATTAAGAATTAGTTTGTTAAGTGATGTTGATCAAATTGTTAGGTTAGCAGGTTTTAAGTCAAATTCTGCCTCCATTGCAAAAAATGGAGTTGCTCAACmTAACATCAGTTGATAATGTTGATGTATCTrTAGACTGGAAT 

973 y 9 e'sv 09 nva'riv'yhsanmt'qms’ievvyvs’ftlvrndgss’ail’dipnfkc'pdtsyclfyk 

3961 TATCAAGAAAGTGTTCAACAGAATGTAGCTCGCATTGTTTATCATTCTGCAAATATGACTCAAATGTCTATTGAAGTTGTGTATGTGAGCrnACGCTTGTGCGTAATGATGGTTCTAGTGCAATACTTGATATACCAAATTTTAAATGTCCAGATACATCTTATTGTTTATTTTATAAG 

1033 PGK'sGVLKFr’GKG TLTSCYD NKHLTWFKVT'cPOFNqPHRL'ATCFVlOQHOVVYPPIKATq 

4141 CCTGGTAAAAGTGGTGTGTTGAAGTTTACAGGGAAAGGGACTTTAACATCATGCTATGATAATAAGAATTTGACCTGGTTTAAGGTTACATGTCCAGATTTTAACCAACCATGGCGTCTTGCTACATGTTTTGTTATTCAACAGCATGATGTTGTTTACCCAGCAATTAAGGCAACACAG 

1093 yenvtfvngppgt’gkttfvy’dtyiskasss'nrfvycapthriv'gdmdekvdgavvvsayn 

4321 TATGAGAATGTCACCTTTGTTATGGGACCACCTGGAACTGGCAAAACAACCTTTGTnATGATACATATCTCTCGAAAGCATCTTCTTCGAACCGTTTTGTCTATTGTGCACCAACACATAGACTTGTTGGAGATATGGATGAGAAGGTTGATGGTGCAGTTGTTGTrTCAGCATArAAC 

1153 ORT’YRNPVWN KOo'sYGVLLc'THN’TLPFIKS AVL’lAOEVSL'lPP HVMiKlL'sKGFKKVVLL 

4501 GATAGGACTTACCGCAATCCAGTTTGGAATAAGGATGATTCTTACGGTGTCCTACTTTGCACACATAACACGTTGCCATTTATTAAAAGTGCTGTGTTGATTGCTGATGAAGTTTCCCTAATACCACCACATGIGATGATTAAAATTTTGTCCATGGGTTTTAAAAAGGTTGTGTTGTTA 

1213 gdp'fqlspvyknh'kvhfkyd'tfy’llqlatq'kry'ltacyrc'pp 9 'ilsafskpycdvgvdlv- 

4681 GGAGATCCATTTCAGnAAGCCCAGTTTATAAAAATCATAAAGTGCATTTTAAATATGATACTTTrrATTTATTACAGCTTGCAACACAAAAACGCTATCTTACAGCATGCTATCGATGTCCACCACAAATCrrAAGTGCTTTTTCAAAGCCTTATTGTGATGTTGGTGTTGATTTAGTT 

1273 SFn'kKPGKFO IIv’sK0LANi'qDf’svLSVLs'kEy'pGYVII.v'nYr'aAVDYAh'0Kg’lGOVTTI 

4861 TCTTTTAATAACAAACCAGGTAAGTTTGATATAATTGTTAGTAAACAGCTTGCCAATATTCAAGATTTTAGTGTACTTTCTGTATTGTCAAAAGAGTACCCTGGTTATGTTATTTTAGTGAACTATAGAGCTGCTGTTGATTATGCAATGCAAAATGGGCTTGGTGATGTGACAACTATT 

1333 dss’ q gttaa n'hil'vlfgasn'fsk’tvnrviv’gcs’rstthiv'vvc'c pelf* h'fop'ilnrpep 

5041 GAnCATCTCAAGGAACAACCGCTGCAAATCATTTGCTGGTATTATTTGGAGCATCCAATTTTTCTAAAACAGTAAATCGTGTTATTGTAGGCTGTTCAAGATCTACTACCCATTTAGTAGTGGTTTGTTGCCCAGAACTTTTTAAGCATTTTCAACCAATTTTAAATTGGCCAGAACCA 

1393 KYR'YFGHEK9's0F’NIlPEVS"sLv’FC0iEFR’HYK'ADPNSKT'RTv'YPG9IAv'vTs’QTL9LYL 

5221 AAGrAFCGFTAnTTGGCAFGGAAAAGCAGTCrGArTTrAATATAATACCAGAAGTGrCATCTCTTGTGTTTTGr6ArAFAGAGmTGGCATTATAAAGCAGACCCCAArTCCAAAACTCGTACAGnTATCCTGGGCAAAnGCTGTTGTAACATCTCAAACT7TACAGTTATACC7T 

1453 gvf’odtgyks'alr'glpkovy’vpp'nrvwkrk'hyp'sye 9 hay"nh 8 Rlfkfii'ott’fgqphfi 

5401 GGAGTCTTTGATGACACTGGCTATAAGTCAGCGCTTCGGGGGCTACCAAAGGATGrTTATGTACCACCTAATTGGGTTTGGATGAGAAAACATTATCCATCTTATGAACAACATGCnATAATATGCAGCGGCTGTTTAAGTTTATAATTGATACAACCTTTGGACAACCATGGTTTATT 

1513 LYSCSKOIKSLKFYVEFDTCYFCSCGEHAICLHRDGRYKCRNCYGGHLISKLVHCKYIOV 

5581 TTGTATAGTTGTAGTAATGATTTAAAGAGmAAAGTTTTATGTGGAATTTGATACTTGTTATTTTrGTAGTTGTGGTGAAATGGCCATTTGCCTAATGCGTGATGGCAATTATAAAIGTAGAAATTGTTATGGTGGAATGTTGATAAGTAAATTAGTTAATTGTAAATATTTAGATGTG 

1573 qke’rvklqda'hoa’icqqfhg'dsh'ealcdav'htk'clylasy'eaa'fkdtihvkykdlcleiq 

5761 caaaaagaaagagttaagttacaagatgcccacgatgctatttgtcagcaatttcatggagattctcatgaagcattgtgtgatgctgtgatgacaaagtgcctgtatttggcaagctatgaagctgcatttaaagatactatacatgtgaaatataaagatctttgcttagaaatacag 

1633 YKl'TSSFVRY OSV HKRYLYR DHG AHYYFRT PRS PMqNVYK YEV GSHAEYS INICTSYEGC 

5941 TATAAAATTACAAGTTCTTTTGTTCGATATGATAGTGTTCATAAACGCTATTTATATCGTGATCATGGTGCCATGTATTATTTTAGAACACCTAGGTCTCCAATGCAAAATGTGTATAAGTATGAAGTTGGTAGTCATGCAGAGTATAGTATAAACATCTGTACTAGTTATGAAGGATGT 

1693 9SF'GKTCTKc’lHl’HClVE9F - HAD’ERFI<EFrL0s’vVKSDYv'E9A’LSPAAKA'LMLTVTKVED 

6121 CAATCATTTGGrAAAACGTGCACAAA£TGTATACATATACATTGCATTGTTGAGCASTTTATGGCA5ATGAAAGATTTAAAGAATTTATACTTGTnCTGTTGTTAA5TCAGATTATGIGGAACAAGCATTGTCACCAGCASCAAAAGCTTTGATGTTAACGSTTACGAAGGTGGAGGAC 

1753 ksf'yiskgvrydl'yoyolsk'svh’rvvnskv'xpl'plysviv'glg’inctvgc'vlpnvphklk 

6301 AAAAGTTTTTATAmCAAATGGTGTTCGCTATGATCTCTATGATTATGACTTGAGTAAATCAGTCATGAGAGTAGTGAACTCCAATGTTAAACCACTCCCTTTGTAnCCGTTATTGTGGGTCTTGGAATTAATTGTACTGTTGGCTGTGTTTTACCTAATGTGCCAATGAAATTAAAA 

1813 dei’litdvpl'stl'rldlqtr’yyi'srptlsh'kns'rrklaga'qvy'dcsrhiyieatgeqpiy 

6481 GATGAATTGCTTATTACAGATGTAGCAnATCAACCTTACGGTTAGATTTGCAAACGTGGTACTATATTTCCTGGCCCACATTAAGTAATAAAAATTCTCGTTGGAAATTAGCTGGAGCACAAGTTTATGATTGTTCTGTGCATATTTATATAGAGGCTACTGGAGAACAACCACmAT 


1873 YLQQGRSESLFELPOT IF S TSRLYlilDHDAAqSFHVSCQlAIETWPSSHHVFSGDFTEVGT 
6661 TATTTACAACAAGGTAAAGGTGAAAGCCTCTnGAACTGCCTGATACACTTTTTTCAACTGGCAGATTGTATAACTTAGACCATGATGCTGCACAGAATTTTAATGTTAAACAGnAGCCATTGAAACAATGCCAAATAATCATCATGTTTTTlIWJGAGATTTTACTGAAGTTGGAACT 
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Fig. 2. Sequence and translation of BEV ORFlb and the 3' end of ORFla. Since ORFlb briefly overlaps with the P gene, the translation of the last block of the 
figure also includes the N-terminal 28 aa residues of the BEV peplomer protein. The nucleotide sequence data in this figure have been assigned the EMBL nucleotide 
sequence database accession number X52374. 


product was precipitated when the C-terminal antiserum was used 
(Fig. 6B). After excision of protein bands from gel and 
scintillation counting, a frameshifting efficiency which ranged 
from 20 to 30% was calculated. 


DISCUSSION 

In this paper we present a substantial part of the polymerase 
sequence of the torovirus prototype, Berne virus. Our data on 
the organization and expression of the 3' part of the BEV 
polymerase gene supplement the previously observed similarities 
between torn- and coronaviruses at the level of gene order and 
translation strategy (4). 

The detected amino acid sequence similarities strongly suggest 
an ancestral relationship between toro~ and coronaviruses, but 
the two virus groups are not very closely related: whereas the 
ORFlb products of relatively distant coronaviruses like the avian 
IBV and the mammalian MHV still contain an average of 56% 
identical amino acids (12), the overall identity between the BEV 
ORFlb product and either of the corresponding coronaviral 
proteins is only about 19%. However, on the basis of the 
sequences included in Fig. 4, the percentage of identical amino 
acids in BEV and EBV/MHV can be calculated to rise to 45 —50% 
for the polymerase and helicase domains (80-90% between IBV 
and MHV) and 40—45% for domain 4 (67% between IBV and 
MHV). 

An analysis of the polymerase domain of toro- and 
coronaviruses does not reveal a closer relationship to either the 
‘Poliovirus-like’ or the ‘Sindbisvirus-like’ superfamily (22). 
However, from a recent survey of the putative helicases of 
positive-stranded RNA plant viruses (29) it may be concluded 
that the toro- and coronaviral domain 3 is more closely related 
to the ‘Sindbisvirus-like’ helicase: almost all conserved amino 
acid residues from the ‘Sindbisvirus-like’ domain (especially from 
the group designated A2; 29) are conserved in the toro- and 
coronaviral domain 3 (Fig. 4), whereas a number of conserved 
residues from the ‘picomavirus-like’ helicase is conspicuously 
absent. 

At the nucleotide level, the sequence of the BEV frameshift 
region illustrates the evolutionary distance between toro- and 
coronaviruses. The sequence in the frameshift region of IBV and 
MHV is well conserved (12). Except for the ‘slippery’ 
heptanucleotide sequence U UUA AAC, no significant sequence 
similarity between BEV and coronaviruses was detected in this 
area. Nevertheless, the RNA structure proposed for both 
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Fig. 3. Proportional dot matrix comparison of the amino acid sequence of the 
predicted BEV ORFlb product and the corresponding polypeptide of IBV (8). 
The plot was generated using the COMPARE program with a window size of 
21 and a stringency of 13. 


coronaviruses seems to be conserved in the BEV frameshift 
region. The predicted pseudoknot structure is essential for 
efficient ribosomal frameshifting in the coronavirus IBV (11) and 
can also be folded in the frameshift-directing sequences of MHV 
(12) and a number of retroviruses (27). The two basepaired 
regions S1 and S2 (28) of the proposed pseudoknot are predicted 
to form a quasi-continuous helix (28; Fig. 5B and 5D). This helix 
is thought to promote frameshifting at the ‘slippery’ sequence 
(11, 30). SI and S2 are connected by loops LI and L2 which 
bridge the deep and shallow grooves of the RNA helix (28; Fig. 
5). In IBV and MHV the L2 region consists of 32 and 33 nt, 
respectively. In BEV two domains which can basepair with a 
part of the loop of the hairpin were identified at 12-16 nt and 
70-74 nt downstream of the hairpin’s base (Fig. 5). Although 
an L2 length of II nt (PKl ; Fig. 5A and 5B) does not pose 


Downloaded from http://nar.oxfordjournals.org/ at University of Exeter on July 22, 2015 




4540 Nucleic Acids Research, Vol. 18, No. 15 


MHV / IBV 


ORF la 

1 2 3 

4 



III 


1 


BEV 

lb 


ORF 1 a 



123 

4 



in 


n 


ORFIb 


.. MHV 595 VIGTTKFYGGWDdMLRrLIkd-4-vI>lGWDYPKCDRAMPlJi-LR-41-yYVKPGGTSSGDATTAfANSVFN-58-SmMILSDDGVVC-41-KGP-HEFCSQHTMLVkmDGDevYLPYPDPSRIL 845 

1) IBV 585 VIGTTKFYGGWDnMLRnLIqg-4-iLMGWDYPKCDRAMPNl-LR-41-lYVKPGGTSSGDATTAyANSVFN-58-SlMILSDDGVVC-41-KGP-HEFCSQHTMLVevDGDpkYLPYPDPSRIL 836 

BEV 515 llGvsK-Y-Glk--fskflkd-9-vf-GsDYt K.CPR tf?-l3fR-37-lln KPGGTSSGDATTAh sNtfyN-51-fl-nf LSDD sfif-34-KGhie EFCS aH-ii-kt aGeY hfLP--srgRlX 744 

*"* *** ** + * * * * ★ * w T TT7T *¥ * 


MHV 9 38 LQSvGaCWCsSQTsLRCGsCIRKPllCCKCaYDHVMsTDHKyVLSvsPYvCnspGC-dvnDVtKLYLGGMSYyCedHKPqySfkL 1022 

2) IBV 929 LQScGvCWCnSQTiLRCGnCIRKPfLCCKCcYDHVHhTDHKnVLSinPYiCsqlGCge-aDViKLYLGGMSYrCgnHKPklSipL 1013 

BEV 845 anfdkvCfcCpnpavsvCeeCyvplpLCayCyYv HVv i snH 8kVedkfkcfCgqdnirelviV--LnnsicmYqCkncvesdrlri 928 


MHV 1218 VQGPPGtGKSHlAIGLAVYyctA-56-TINALPElvtDIivVDEVSML-16-yVYiGDPAQLP-30“FLgtCYRCPKEIVDTVS-72-LGLqtQTVDSaQG-19-vNRFNVAiTRAKkGILc 1512 
IBV 1210 VQGPPGsGKSHfAIGLAVYrS8A-56-TINALPEv8cDIllVDEVSML-16-vVYvGDPAQLP-29-FLakCYRCPKEIVDTVS-75-IXiLnv<5TVDSsQG-19-iN _ RPNVAlTRAKrGILv 1506 

BEV 1099 VmGPPGtGKttfvyd--tYlskA-53-ThNtLPfiksavliaDEVSli-15-vVHGDPfQLs-25-yLtaCYRCPpqIlsafS-64-LGdvt-TiDSsQG-19-vNRviVgcsRsthlvvv 1374 
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MHV 2322 HwYGsfnq- -kiiGGLHllIGl-38-TViDLLLDDFvdivk- 7-SKWnVnvDf kdfqFMlVcn-eekvmTfYP-29-lwNYG-kpITLP-5-NVAKYTQLCQYLStTTlaVPaNMRVlHlGA 
4) IBV 2274 HilYGevdk--pqlGGLHtvIGm-38-TVvDLLLDDFlellr-10-SKWtVsiDyhsinFMtWfe-dgsikTcYP-27-ipNYG-vgITLP-5-NVAKYTQLCQYLSkTTicVPhNMRVmHfGA 

BEV 1921 Hvf sGdf tevgtd iGGvHhwal -29- TlvDv canqlyekvk- 8- SKVi fVn iDf qdvqFMvfangeddiqTfYP-23-l kNYG qnptfmP-3-Nf AKYTQiC t fi -qdhvkVarNalVwHlGA 


GSDKGVAPGSaVL-39-DLi.ISDMY-28-LALGGSvAiKi.TEf S-17-FCTn-VNASSSEgFLI 2641 
GSDKGVAPGStVL-39-DLvISDMY-32-LALGGS£AvKvTEtS-17-FCTa-VNASSSEaFLI 2598 
agvdGcaPGdiVL-32- nLivSDiY -16- LALGGt. ivfKtTEsS-17-FfTagVNt SSSE v Fw 2209 


Fig. 4. Positions and sequences of the conserved domains in the ORFIb product of the POL genes of BEV, IBV Beaudette (M42; 8) and MHV-A59 (12). Amino 
acids identical in MHV and IBV are shown in capitals. Amino acids in BEV which are identical in BEV, MHV and IBV are also shown in capitals. Amino acids 
which are identical or conserved in all three viruses are underlined in the BEV sequence. Domain 1: amino acids described to be conserved in the majority of positive- 
stranded RNA viruses (22) are indicated by asterisks. Domain 2: conserved cysteine and histidine residues are marked by asterisks. Domain 3: amino acids conserved 
in ‘Sindbisvirus-like’ RNA plant viruses (group A2; 29) are indicated by asterisks. 
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Fig. 5. Predicted secondary and tertiary RNA structure of the BEV ORFla/ORFlb overlap region. A) Predicted structure using the most upstream region (nt 1099-1103) 
which is complementary to the loop of the hairpin to form an RNA pseudoknot (PK1). The ‘slippery’ sequence U UUA AAC is indicated by a dashed line. The 
ORF la termination codon is underlined. B) Schematic representation of PK1. The ‘slippery’ sequence U UUA AAC is boxed. The basepaired stem structures SI 
& S2 and the connecting loops LI & L2 (28) are indicated. C) Alternative model using the sequence from nt 1157-1161 for pseudoknot folding (PK2). The Ddel 
site used to generate construct pBSFS<5 is indicated. D) Schematic representation of PK2. The 49 nt box indicates the possible internal stem-loop structure in L2 
as shown in panel C. 


any sterical constraints upon pseudoknot formation (28), the 
involvement of the more downstream complementary region 
(PK2; Fig. 5C and 5D) appeared to be an interesting alternative: 
the resulting 69 nt L2 contains possibilities for internal folding 


(Fig. 5C), which could play a role in the formation or stabilization 
of the pseudoknot structure. 

An attempt was made to establish whether the PK1 or PK2 
is involved in the formation of the BEV pseudoknot structure. 
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Fig. 6. Analysis of BEV ribosomal frameshifting in vitro and in vivo using construct pBSFS. A) Schematic representation of the pBSFS construct. The black box 
represents the T7 RNA polymerase promoter in pBS. Open boxes indicate BEV sequences; N- and C-terminal MHV M sequences are represented by differently 
hatched boxes. The composition of the two predicted translation products is shown. B) SDS-PAGE of pBSFS translation products immunoprccipitatcd by the N- 
terminal (aN) or C-terminal (aC) MHV M antisera. In the lanes marked * + ’ in the in vitro panel, pBSFS transcripts were added to a rabbit reticulocyte lysate. 
The lanes marked ' —' show a control where no RNA was added to the translation system. The in vivo panel shows the results of an experiment with a pBSFS- 
transfccted, vTF7-3-infected Hela cell lysate (lanes marked '+') and an untransfected vTF7-3-infected control lysate (lanes marked * — ’). 


For this purpose construct pBSFSS was generated (using a Ddel 
site at position 1145-1149 in the sequence; Fig. 5C) which 
contains a smaller part of BEV ORFlb and thereby lacks the 
downstream complementary domain of PK2 (nt position 
1155-1161 in Fig. 2). However, the deletion did not result in 
a reduction of the frameshifting efficiency (data not shown), 
suggesting that PK2 is not required for efficient frameshifting. 
It should be noticed, however, that our present data do not 
formally exclude a role for this structure in translational 
frameshifting. If both pseudoknots are involved in the process, 
one structure could maintain the frameshifting efficiency in the 
absence of the second pseudoknot. 

In conclusion, our data on the expression, organization and 
translation product of the 3' part of the BEV polymerase gene 
provide additional evidence that toro- and coronaviruses are 
evolutionarily more closely related to each other than to any other 
family of positive stranded RNA viruses. If the toroviruses are 
to be classified as a separate virus family, they must be considered 
members of a new superfamily of positive-stranded RNA viruses 
(6, 7), together with the coronaviruses. Alternatively, our data 
would justify a taxonomic proposal in which corona- and 
toroviruses form two separate genera of a new family of positive- 
stranded RNA viruses. 
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